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Abstract 
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1 Introduction 


Schwinger’s original nonet mass formula Jl|] (here the symbol for the meson stands 
either for its mass or mass squared), 

(4 K — 3rj — tt)(3 77 ' + n - 4 K) = 8(K - tt) 2 , (1) 


and the Sakurai mass-mixing angle formula 0, 


tan 2 d 


4 K — 3r/ — it 
3rj' + 7 r — 4 K ’ 


( 2 ) 


both relate the masses of the isovector (77), isodoublet (K) and isoscalar mostly octet 
( 77 ) and mostly singlet (77 ) states of a meson nonet, and the nonet mixing angle (d). 
We alert the reader that, although we use notation suggestive of masses below, each 
formula is to be reprised in terms of mass or mass squared values in this introductory 
section. 

The relations (1) and (2) are usually derived in the following way: For a meson 
nonet, the isoscalar octet-singlet mass matrix, 


M = 


M$s Mg 9 
Mgy M99 


(3) 


is diagonalized by the masses of the physical 77 and 77 states: 

M = 


cos d sin 9 
— sin 6 cos 9 


7 ] 0 

0 7/ 


cos 9 
sin 9 


■ sin 9 
cos 9 


cos 2 9 77 + sin 2 9 r\ sin 9 cos 9 (rj' — rf) 
sin 9 cos 9 (77 — 77 ) sin 2 9 77 + cos 2 9 77 


(4) 


where 9 is the nonet mixing angle, which is determined by comparing the correspond¬ 
ing quadrants of the matrices (3) and (4), by any of the three following relations: 


tan 2 9 = 


tan 2 9 = 


Mi 


88 


sin 29 = 


r( — Mg 8 ’ 

rj ~ M 99 

M 99 - 77 ’ 

2Mgg 


(5) 


( 6 ) 


(7) 


77 — 77 

ft is easily seen that Eqs. (5) and ( 6 ) are identical, since, due to the trace invariance 
of M, 

77 + 77 ' = M 88 + M 99 , ( 8 ) 

and therefore, M 88 — 77 = q — M 99 , and q — M 88 = Mg 9 — 77 . Eliminating 9 from 
(5),(7), or ( 6 ),(7), with the help of sin 2d = 2 tand/(l + tan 2 d), leads, respectively, to 


(77 - M 88 )(M 88 - 77 ') = Ml g, 


(9) 
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(M 99 - ri)(rf - M m ) = M|g, 


( 10 ) 


which again are identical, through (8). 

We note that, under the quark model inspired conditions 


M 8 8 — 


4 K 


it 


M 99 — 


2 K + it 


Ms 9 — 


2\/2 


{K -tt), 


( 11 ) 


where the first of the three relations in (11) is the standard Gell-Mann-Okubo mass 
formula ||, Eqs. (9),(10) have the effectively unique solution 


n 

V 


2 K — it, ij = it, 6 
7T, 77 = 2 K — it, 6 


1 

arctan —= = 

V2 

— arctan \/2 = 


35.3°, or 


-54.7°, 


( 12 ) 

(13) 


which corresponds to the “ideal” structure of a meson nonet. 

For all well established meson nonets, except the pseudoscalar (and, we expect, 
scalar) one(s), both linear and quadratic versions of Eqs. (5)-(10) are in good agree¬ 
ment with experiment. For example, for vector mesons, if one ass um es the validity of 
the Gell-Mann-Okubo formula u / 8 = (4/1 * — p)/3, then one obtains from Eq. (9) with 
the measured meson masses fj, M 89 = —0.209 ± 0.001 GeV 2 in the quadratic case, 
and —0.113 ± 0.001 GeV in the linear case. Note that this is entirely consistent with 
—0.196 ± 0.005 GeV 2 and —0.118 ± 0.003 GeV, respectively, which follow from the 
third element of (16). For tensor mesons, a similar comparison gives —0.305 ± 0.020 
GeV 2 vs. -0.287 ± 0.015 GeV 2 , and -0.105 ± 0.008 GeV vs. -0.104 ± 0.005 GeV, 
respectively. 

However, for the pseudoscalar nonet, one obtains from Eq. (5) with meson masses 
squared, 6 ss — 11 °, in sharp disagreement with experiment, which favors the 77-/7 
mixing angle in the vicinity of — 20 ° 1 i 1 - Although using linear meson masses 
in Eq. (5) does give 6 sa —24°, in better agreement with data than its mass-squared 
counterpart, the value of M 89 , as given by (9), is now —0.165 ± 0.004 GeV, vs. 
—2^/2/3 (K — it) = —0.338 ± 0.004 GeV. This emphasizes that neither Schwinger’s 
nonet mass formula nor the mass-mixing angle relations (including Sakurai’s) (5)-(7) 
hold for the pseudoscalar nonet. 

ft is well known, however, that the pseudoscalar (and, probably, scalar) mass 
spectrum does not follow the “ideal” structure, Eq. (11), since the mass of the 
pseudoscalar isoscalar singlet state is shifted up from its “ideal” value of (2 K + 7t)/3, 
presumably by the instanton-induced’t Hooft interaction |7| which breaks axial U(l) 
symmetry [H Pi |l0i . However, the use of M 99 = (2/1 +7t)/3 + A, 4^0, in Eqs. (5)-( 8 ) 
will again lead to Schwinger’s formula (9), which does not hold for the pseudoscalar 
mesons, as just demonstrated. [In fact, the structure of this formula does not depend 
at all on M 99 , as seen in (9).] Therefore, instanton, as well as any other effects which 
may shift the mass of the pseudoscalar isoscalar singlet state, cannot constitute the 
explanation of the failure of Schwinger’s quartic mass and the Sakurai mass-mixing 
angle formulae for the pseudoscalar nonet. We believe, however, that the following 
analysis can resolve this problem. 
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2 Pseudoscalar meson mass squared matrix 

Here we suggest the following form of the mass squared matrix for the pseudoscalar 
mesons: 


p M 2 _l( ifl-K 2 - /> 2 -2 V2 (fl-K 2 - jy‘) \ 

1 3 f -2V2 UkK 2 - /> 2 ) 2/J/i' 2 + /> 2 + 3/ 2 A J 


(14) 


where /’s are the pseudoscalar decay constants defined below, and A stands for 
the sum of all possible contributions to the shift of the isoscalar singlet mass (from 
instanton effects, l/lV c -expansion diagrams, gluon annihilation diagrams, etc.). 

Indeed, the form of such a mass squared matrix is determined by the form of 
Gell-Mann-Okubo type relations among the masses of the isovector, isodoublet, and 
isoscalar octet and singlet states (which in our case are Eqs. (17)-(19) below), since 
this matrix must be equivalent to that of the form (3), which in the case we are 
considering is 


f 2 M 2 


fivis fsf9vi 9 \ 
fM fWm )' 


(15) 


and which we assume is diagonal in the basis of the physical r/ and r/ meson masses 
and decay constants: 

PM 2 = ( f ^ ) . (15) 

The equivalence of the matrices (14) and (16) is guaranteed by the validity of the 
following relations: 


fsVss = m u (uu) + m d (dd) + 4 m s (ss) 




(17) 


/sMsq = --^r \ m n{uu) + m d (dd) - 2m s (ss ) ] N 


3 L 


fldm = fgA ~ o m n(uu) + m d (dd) + m s (ss) = f$A + 


{/y - s 2 k k 2 ) , 
2 flu 2 + py 


with 


K = 


(Jl ±) 2 + (K°f 


(18) 

(19) 

( 20 ) 


as suggested by Dmitrasinovic [JT0[] , on the basis of the (precise) Gell-Mann-Oakes- 


Renner formulae which relate the pseudoscalar masses and decay constants to the 
quark masses and condensates [11]: 


fy 2 = - 

m u (uu) + m d (dd) 

/iM 2 = - 

m u (uu) + m s (ss) 

/i(A-") 2 = - 

m d (dd) + m s (ss) 


( 21 ) 

( 22 ) 

(23) 
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(We ignore Dashen’s theorem violating effects [|T2|, and only approximately take into 
account isospin violating effects via (20), as we are not concerned here with accuracies 
better than 1%.) 

Note that in the limit of exact nonet symmetry, 


fn — fK — fs — /j) = /) 


(24) 


the mass squared matrix (14) reduces to (3) with the matrix elements defined in (11) 
and an additional contribution in the isoscalar singlet channel (A). 


3 Modified Gell-Mann—Okubo mass formula 

Here, we shall only explicitly demonstrate the validity of the modified Gell-Mann- 
Okubo formula (17) (the remaining relations (18) and (19) may be checked in a similar 
way). This formula may be obtained in standard chiral perturbation theory [[13|], as 
follows. 

Standard chiral perturbation theory leads to the following expressions for the 
pseudoscalar meson masses and decay constants ||T3|: 


7T 2 = 


K 2 = 


2 m n B ^1 + /j, n - -p 8 + 2 m n K 3 + (2 m n + m s )K^j , 

K + m s)B ^1 + + (m n + m s )K 3 + (2m n + m s )K^j , 

2 / 4 2 \ 

hss = g( m n + 2m s )B (1 + 2fi K — -/i 8 + -{m n + 2m s )K 3 + {2m n + m s )KA 

+ 2 m n B i-Hn + -hk + -hs ) + B(m s - m n ) 1 2 K 5 , 
fir f 2 hk 2m„ /\ (i T (nr n T 2m s )K 7 ^ , 

Ik = f (1 - -/A- - -Pa' - -Ps + (m« + m s )AT 6 + (m n + 2m s )K 7 j , 


/ ( 1 - 3/ift- + -(m n + 2m s )K 6 + (m n + 2m s )K 7 ) , 


(25) 


where p’s are chiral logarithms, and the constants A'; arcQ proper combinations of 
the low energy coupling constants L lr It then follows from these relations that the 


standard Gcll-Mann-Okubo formula is broken in first nonleading order [13 


A 


GMO 


= 4A' 2 — 3//g g — 7r 2 


= 4 B 


m n (ii k + p 8 - 2 hr) + 2m s (p 8 - fj, K ) + ■■■ , 


(26) 


where ... indicates higher order terms. 


1 We have extracted factors absorbed into K 4 , K 5 and K 7 in ref. |L3|, to make uniformly explicit 

the quark mass factors and the overall scale. 
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However, the modified Gcll-Mann-Okubo formula remains valid in this order, and 
is violated only by second (non-leading) order SU(3)-flavor breaking effects: 


A gmo = 4 ( 4 /i A ' 2 3 - 3/|^I 8 - fW) 


f 

= AB (m s - m n ) 




- (m s - m n ) ( |/i 3 + -Ji 5 + -K -; 


+ 


(27) 


since the second factor in the first bracketed term on the r.h.s. of (27) must vanish 
in the SU(3)-flavor limit. 

This analysis shows, in an essentially model independent way, that the modified 
Gell-Mann-Okubo formula (17) must be a more accurate relation among the octet 
pseudoscalar mesons than the standard one,Q and therefore, use of the form of the 
mass squared matrix (14) is fully justified.^ 


4 The Schwinger and Sakurai formulae reexam¬ 
ined 


Starting with the mass squared matrix (14), the considerations which lead to Eqs. 
(5)-(10) above, will now lead, through (17)-(19), to the following two relations, 


tan 2 6 


*f 2 K K 2 - 3/,V - /> 2 
3/^' 2 + /> 2 - 4 /^’ 


(28) 


(4 - 3 /,y 2 - /„v) ( 3 /yy + ;y 


4 f 2 K K 2 


) = 8 (f%K 2 - f 2 y ) 2 , (29) 


where 9 is the mixing angle for the valence quark wavefunctions which produces the 
mass eigenstates from the octet-singlet basis, and which we refer to as “the Schwinger 
nonet mass, and Sakurai mass-mixing angle (respectively), formulae reexamined”. 

In contrast to f n and fx, the values of which are well established experimentally 

a. 


V2f„ = 159.8 ± 1.6 MeV, = 130.7 ± 0.3 MeV, ij-= 1.22 ±0.01 (30) 

Jn 

the values of f v , f > and 6 are known rather poorly. We now wish to calculate 
these values using the relations (28),(29), and compare the results with available 
experimental data. It is obvious that the two relations are not enough for determining 

2 A search for a relation of a more general form, Af^K 2 = 3fgr]g 8 + /“7r 2 , which would hold in 
the first nonleading order of standard chiral perturbation theory, results in a = 2. 

3 It has been suggested in the literature that the pseudoscalar decay constants should enter rela¬ 
tions like (17)-(19) in the first rather than second power |IJ|. As discussed above, such relations are 
expected to be less accurate than ours, according to chiral perturbation theory. 
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the three unknowns. However, the additional relation, independent of (28),(29) (the 
trace condition for (14),(16)), 

fW + f}v 2 = 2 f 2 K K 2 ± /g A (31) 

introduces yet an additional unknown, A. We therefore develop another independent 
relation among f n , f > and 9, as follows. 

The light neutral pseudoscalar decay constants are defined by the matrix elements 

<0| E<A(°bA 5 A(0)l Pip)) = ti jP fpp^ (32) 

q 

where ij} q , q = (u, d, s ) are the quark held operators, and we define the wave function 
of a neutral pseudoscalar meson P in terms of the quark basis states qq as: 

! P ) = E C JEA q = u,d, s, (33) 

q 


where for P = it 0 , c® = l/\/2 = — c d , c® = 0, for P = r/ 88 , c® = c® = 1/vA 
c® = —2/\/6 and for P = 7799 , c^ = c 9 d = c 9 s = l/\/3. 

The pseudoscalar decay constants dehned in (32) can now be expressed as 

fp = E( c f) 2 /^> (3 4 ) 


where we have introduced the auxiliary decay constants f qq dehned as the decay 
constants of the qq pseudoscalar bound states having the masses M(qq). In the isospin 
limit, fuu = fdd — fud — f-n- 0 — fn+ ■ Using this approximation, and evaluating the 
appropriate matrix elements leads to the following relations: 


fv 

u 


cos 


9 — y/2 sin 9 ' 


Vz 

' sin 9 + y/2 cos 9 ' 


fn + 


fn + 


sin 9 ± v^cos^ 

2 

) faa, 

(35) 

V3 , 

cos 9 — y/2 sin 9 N 

to 

to! 

(36) 

Vs > 


Now we have four equations, (28),(29),(35),(36), which allow us to determine the 
three unknowns, f r] , / / , 9, as well as the additional quantity introduced, namely, f S g. 
The solution to these four equations is 


yi- = 1.085 ±0.025, 

Jn 

f, 

yL = 1.195 ±0.035, 

In 
f - 

= 1.280 ±0.060, 

In 

g = (-21.4 ±1.0)°. 


(37) 

(38) 

(39) 

(40) 


[The 7 T and K electromagnetic mass differences, and the uncertainties in the values 
of f n and fx, (see (30)), are taken as a measure of the uncertainties of the results.] 
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Before comparing the solution obtained with experiment, let us also calculate the 
values of fa and fa which are obtained from (35),(36) in the no-mixing case (9 = 0): 


1 2 

J 8 3 JtI- T g/sS) 

(41) 

2 1 

fa r^fa T g/sS' 

(42) 

Therefore, as follows from (39),(41),(42), 


fo 

y- = 1.185 ±0.040, 
fa 

(43) 

j- = 1.095 ±0.020. 

/ 7T 

(44) 


The rj-q' mixing angle, as given in (40), is in agreement with most of experimental 
data which concentrate around — 20 ° [§ s 1- Also, the values for fa/fa, fa/fa and 
9 are consistent with those suggested in the literature, as we show in Table I. 


Ref. 

fa/fa 

fa/fa 

0, deg. 

This work 

1.185 ±0.040 

1.095 ±0.020 

-21.4 ± 1.0 

[ 6 ] 

1.11 ±0.06 

1.10 ± 0.02 

-16.4 ± 1.2 

[15,16] 

1.25 

1.04 ±0.04 

-23 ±3 

[17] 

1.33 ±0.02 

1.05 ±0.04 

-22 ±3 

[18] 

1.12 ±0.14 

1.04 ±0.08 

-18.9 ±2.0 

119] 

1.38 ±0.22 

1.06 ±0.03 

-22.0 ±3.3 

[ 20 ] 

1.254 

1.127 

-19.3 


Table I. Comparison of the values for fa/ fa, fa/fa and 9, calculated in the paper, 
with the results of the papers referenced. 


5 Comparison with data 

We now wish to compare the values obtained above for the ratios fa/fa, fa/fa, and 
for the ij-r] mixing angle with available experimental data. We shall first consider in 
more detail the well-known 77 , 77 —> 77 decays, for which experimental data are more 
complete than those for other processes involving light neutral pseudoscalar mesons, 
and then briefly mention the 7 , 7] —>• 7 r + 7 r“ 7 , and J/t[j —> 7 / 7 , 7] 7 decays. 


5.1 


7/, 7/ 


77 decays 


For these processes, the inclusion of both the 77-77 mixing and the renormalization of 
the octet-singlet couplings, which leads to the predicted amplitudes 


^77(0) — 


aN r 






3V0nfa 

CtN r . 


'fa n 

— cos 9 


2y/2^j- sin 9 
fa t 


fa 




fa 


(45) 


















The values of these amplitudes, as extracted from data on widths, are [p 


iv, 77 (0) = 0.024 ± 0.001 GeV 
F Vry (0) = 0.031 ±0.001 GeV 


Calculation with the help of Eqs. (30),(40),(43)-(46) yields 

±^ 77 (0) = 0.025 ± 0.001 GeV -1 , 
i± 77 (0) = 0.030 ±0.001 GeV" 1 , 

in excellent agreement with (47). 


(47) 


(48) 


5.2 r] —> 7T + 7r 7 decays 

These, as well as the P° —* 77 , processes were extensively studied by Venugopal and 


Holstein |19| in chiral perturbation theory. The analysis of experimental data for 


both of these processes done in ref. [19] yields 


j- = 1.38 ±0.22, ^ = 1.06 ±0.03, 9 = (-22.0 ± 3.3)°, 

J 7V J 7T 


in good agreement with our Eqs. (40),(43),(44). 

5.3 J/-0 —> r/ 7 , 7/7 decays 


These processes were studied by Kisselev and Petrov |T8|. The values of /§/ 
fg/f n and 9 extracted in ref. ]18| from the experimentally measured P° —>■ 77 and 
J/tj} —> 777,77 7 widths, as given in the last column of Table I of ref. JT8[, which 
corresponds to conventional mass-mixing angle relations, are 

j- = 1.12 ±0.14, ^ = 1.04 ±0.08, 9 = (-18.9 ±2.0)°, 

J 7r Jn 

again in good agreement with our Eqs. (40),(43),(44). 

Thus, the three values of fs/fn, fv /and 9 agree with experiment (at least, as 
far as the processes considered above are concerned). 


5.4 Comparison with CELLO and TPC /27 results 

As to the remaining /,,//, r , f> /f n ratios also calculated in the paper, the experi¬ 


mental values of them, as extracted from data by the CELLO |2T[ and TPC /27 
collaborations, are, respectively, 

fv 


^ = 1.12 ± 0 . 12 , 
J7T 


k 

u 


= 1.06 ± 0 . 10 , 


(49) 

(50) 
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and 


u 

u 

k 

u 


1.09 ±0.10, 
0.93 ±0.09. 


(51) 

(52) 


While the value calculated for f v / f n , Eq. (37), clearly agrees with both experimental 
values (49) and (51), the value calculated for f v r /f n , Eq. (38), only marginally agrees 
with (50), and disagrees with (52) by almost 3 standard deviations. 

To clarify the difference from the results referred to above, we recall that the 
values of f v /f n and f v '/fn are extracted by both CELLO and TPC/2y from exper¬ 
imental data on transition form-factors TW/)(0, — Q 2 ), assuming that the pole mass 
A??!??'), which parametrizes their fits to the data, can be identified with 2n^2 f^y 
Then, these pole fits to the data are presumed to join smoothly, as Q 2 —> oo, to the 
perturbative QCD predictions for T^q(0, — Q 2 ) |23|, i.e., these fits would then agree 
with both the QCD asymptotic form ~ 1 /Q 2 and its coefficient. 

However, the values of both and f > quoted by the two groups, (in MeV) 
(94.0 ± 7.1, 89.1 ± 4.9) 0, and (91.2 ±5.7, 77.8±4.9) 0, respectively, are all close 
to M(p)/(27rv / 2) ~ 86.5 MeV, thus indicating a possible connection with the vector 
meson dominance interpretation of A ,is ps M(p) in the range of Q 2 investigated, 
which could obviate the assumptions referred to above. 

On the other hand, as remarked by Klabucar and Kekez |24|, their own bound- 


state calculation, as well as the model independent calculation by Gasser and Leutwyler 
0 , testing the Goldberger-Treiman relations by Scadron [0, and the calculation 
by Burden et al. 0] , all agree that both f v and f> should be noticeably larger than 
f n . Our results /, 7 1-1 fir, f v 1.2 f w are in agreement with this. 

It therefore seems possible that the discrepancy with other determinations may 
be due an inability to extract the values of f v and / / from the transition form- 
factors TWO, — Q 2 ) Sufficiently accurately, at least in the range of Q 2 investigated 
so far. That this may indeed be the case is indicated by the experimental value 
/> = 84.1 ± 2.8 MeV 0 , extracted by the same method, which is again close 
to M(p)/27ry / 2. The central value of this f n o, 84.1 MeV, is ~ 10% below the well 
established value given in Eq. (30), 92.4±0.2 MeV. Such a large discrepancy cannot be 
explained by, e.g., small isospin violation, indicating the possibility that the extracted 
values for both f v and / / may well have been underestimated also. 

6 Concluding remarks 

As a lagniappe, we note that Eq. (31) may be combined with our extracted values 
for /g, f v , f r > to obtain the value of A : 

A = 0.78 ± 0.12 GeV 2 . 

This is consistent with the usual value 0.73 GeV 2 determined by the trace condition 


for (14),(16) without /’s |^7 
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Let us briefly summarize the findings of this work: 

i) We have found that chiral perturbation theory suggests that the natural object 
to study is the pseudoscalar mass squared matrix modified by the inclusion of the 
squared factors of the pseudoscalar decay constants, rather than the mass squared 
matrix itself. 

ii) We have shown that this modified mass squared matrix leads to new Schwinger’s 
quartic mass and Sakurai mass-mixing angle relations for the pseudoscalar meson 
nonet. 

iii) We have used these new relations to calculate the pseudoscalar decay con¬ 
stants and mixing angle. We have demonstrated that, except where questions may 
be raised regarding the reliability of the extraction of the relevant quantities from di¬ 
rect experimental data, the results obtained are in excellent agreement with available 
data. 

We thank Philip R. Page for useful conversations. 
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